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Y^ ; Abstract 

^N I We explore the potential offered by large-mass Liquid Argon detectors for deter- 

C^^ ' mination of the sign of Am|]^, or the neutrino mass hierarchy, through interactions of 

atmospheric neutrinos. We give results for a 100 kT sized magnetized detector which 

^ . provides separate sensitivity to i/^ , i/^ and, over a limited energy range, to u^ , I'e ■ 

We also discuss the sensitivity for the unmagnetized version of such a detector. After 
including the effect of smearing in neutrino energy and direction and incorporating the 
relevant statistical, theoretical and systematic errors, we perform a binned x^ analysis 

5^ I of simulated data. The x^ is marginalized over the presently allowed ranges of neu- 

trino parameters and determined as a function of ^13. We find that such a detector 
offers superior capabilities for hierarchy resolution, allowing a > 4o" determination for 
a 100 kT detector over a 10 year running period for values of sin^ 2^13 > 0.05. For an 
unmagnetized detector, a 2.5a hierarchy sensitivity is possible for sin^ 2^13 = 0.04. 
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1 Introduction 

Parameters related to neutrino masses and mixings play a fundamental role in efforts to 
construct a viable theory beyond the Standard Model. Over the last decade, our knowledge 
of these parameters has increased at an unprecedented pace, due to crucial results from solar, 
atmospheric, reactor and accelerator based neutrino oscillation experiments. (For a recent 
review, we refer the reader to [1] and references therein.) 

Among the important but as yet unanswered questions, one of the foremost is the determi- 
nation of sign(Am3]^) or the hierarchy of neutrino massed, presently unconstrained by available 
data. If sign(Am|^) > 0, then we have the mass pattern, ma ^ m2 ^ mi, which is similar 
to that of the charged leptons. This is called the normal hierarchy (NH). If sign(Am|^) < 0, 
then the mass pattern is m2 > mi ^ ma. This is called the inverted hierarchy (IH). In recent 
papers, the prospects for progress towards a resolution of this question offered by atmospheric 
neutrinos have been explored [2-10]. With most of their flux below 10 GeV, these neutrinos 
traverse distances up to ~ 12000 km in the earth's matter enroute to a detector. This exposes 
them to appreciable resonant matter effects which occur (primarily) for energies between 2—10 
GeV and distances between 4000—12500 km inside the earth [5,11,12]. Until very long baseline 
experiments using /3-beams or neutrino factories [13] are built, atmospheric neutrinos permit 
us to exploit these effects, albeit in a slower and less spectacular fashion which calls for a care- 
ful analysis of accumulated effects over many bins in energy and angle. Such an approach was 
recently pursued in [6,10] for megaton water Cerenkov detectors {e.g. HK [14,15], UNO [16] or 
MEMPHYS [17]) and magnetized iron detectors {e.g. I NO [18]). The salient features emerging 
out of the above analyses are 

• Both muon and electron events arising from atmospheric neutrinos have sensitivity to 
sign of Am|^ for sufficiently large values of sin^ 2^13 . 

• One of the major factors responsible for a reduced sensitivity is the finite energy and 
angular smearing of a detector. 

• The muon events are more sensitive to smearing than electron events. 

• Magnetized iron calorimeter detectors have the advantage of having charge identification 
capability but they are sensitive only to muon events 

• Water Cerenkov detectors have the advantage that they are sensitive to both muon and 
electron events. On the other hand they do not have charge sensitivity. 

The above features indicate that a detector which is sensitive both to muon and electron 
events as well as their charges will be ideal for probing hierarchy of neutrino masses. An 
important class of future detectors using Liquid Argon as their active medium, may provide 
this kind of set up specially if there can be magnetized versions of these. In this paper we 
study the hierarchy sensitivity of Liquid Argon detectors for cases with and without charge 
sensitivity. 



"'^We use the convention Am?; = m^ — m?. 



2 Basic Characteristics of Liquid Argon Detectors 

Liquid Argon detectors are time projection chambers with fine-grained tracking and total 
absorption calorimetry. Ionization electrons resulting from the passage of an energetic particle 
through the medium are detected by drifting their paths over several meters to wire planes. 
The orientation of these planes is designed to reconstruct the time, length and position of each 
path by recording multiple snapshots of the electrons, from which a bubble-chamber like image 
is constructed. The viability of this technology for a small detector has been convincingly 
demonstrated by [19], and intensive efforts are underway for the upgraded development of 
large mass detectors [20-22] to fully exploit the promise of this technique. 

Essentially, the technique allows us to detect charged particles with good resolution over 
the range of MeV to multi GeV. It uses well studied dE/dx measurements in the medium to 
separate electrons, muons, pions, kaons and protons from each other. It is also possible to 
separate the light from electrons vs that from neutral pions with high efficiency. Magnetization 
over a 100 kT volume has been deemed possible [23,24]. As mentioned earlier, we present 
results for both the magnetized and unmagnetized case. 

We assume the following energy resolutions over the ranges that will be relevant to our 
calculations [21]: 
For the GeV electrons that we will be interested in, 



ajE, = 3%/V^e (1) 

where Eg^ is the electron energy and o"e is the energy resolution of the electron. All values of 
energy are expressed in GeV. For the muon neutrino we take the energy resolution in terms 
of the muon energy to be 

ajE^ = 15% (2) 

Additionally, for hadronic showers , which account for a significant fraction of the uncer- 
tainty in the determination of the primary neutrino energy, we assume 

(^hadlEhad = 30%/ VEhad (3) 

where E^ad is the hadron energy and ahad is the energy resolution of the hadron. 

Our procedure for inferring the neutrino energy from the measured charged lepton and 
hadronic energies incorporates a knowledge of the average rapidity for charged current cross- 
sections in this energy region. The energy resolution in terms of the neutrino energy is related 
to the leptonic and hadronic energy resolutions as follows: 



a^/E^ = V (1 - yf{(yieplEiepf + y^{ahad/Ehady (4) 

The rapidity is defined as y = Ehad/E,^, where E^, = Eiep + Ehad is the energy of the neutrino. 
Therefore, the energy resolution in terms of the neutrino energy is given by 



^ujE,^ = v/(l-l/)(0.03)Vi?.,+y(0.3)VE,, (5) 



and 

<^ujE^^ = y/(l-|/)2(0.15)2 + y(0.3)V^.^ (6) 

for electron and muon neutrinos respectively. 

In our computation, we take the average rapidity in the energy region of our interest 
(i.e. in the GeV range) to be 0.45 for neutrinos and 0.3 for antineutrinos [25]. The angular 
resolution of the detector is taken to be ag = 10°. The energy threshold and ranges in which 
charge identification is feasible are Ethreshoid = 800 MeV for charge identification of muons 
with high efficiency and Eeiectron = 1 — 5 GeV for charge identification of electrons with high 
efficiency. Charged lepton detection and separation (e vs n) without charge identification 
is possible for Eiepton > few MeV. Also, for electron events, it is expected to have a 20% 
probability of ~ 100% charge identification in the energy range 1-5 GeV. 

The next section provides a brief description of our calculational and numerical procedure, 
prior to our discussion of the results. 

3 Calculational and Numerical Procedure 

We follow the procedure described in [10], and refer the reader to the description provided 
there for any details that may be omitted here. 

Our calculation of the atmospheric electron and muon event rates uses the neutrino oscil- 
lation probabilities corresponding to the disappearance channels P^^ and Pee and appearance 
channels P^e and Pe^ for both neutrinos and antineutrino^l- We have modeled the density 
profile of the earth by the Preliminary Reference Earth Model (PREM) [26] in order to nu- 
merically solve the full three fiavour neutrino propagation equation through the earth. The 
total CC cross section used here is the sum of quasi-elastic, single meson production and deep 
inelastic cross sections [27-30]. For the incident atmospheric neutrino fiuxes we use the tables 
from [31]. 

For our analysis in Liquid Argon detectors, we consider an exposure of 1 Mt yr (100 kT 
X 10 years). We look at the neutrino energy range of 1 - 10 GeV and the cosine of the zenith 
angle (6) range of -1.0 to -0.1. These ranges are divided into 9 bins in energy and 18 bins in 
zenith angle. The /i~ event rate in a specific energy bin with width dE and the solid angle 
bin with width dQ is expressed as : 



d^N^ _ 1 



dl^dE 2tt 



d^^^ ]p +( d^^- IP 
dcosMEJ '^'^^ VdcosME^ '^'^ 



CTCC Deflf (7) 



Here <l>^,e are the atmospheric fiuxes (z/^ and Ue), o"cc is the total muon-nucleon charged 
current cross-section and D^s is the detector efficiency. The /i"*" event rate is similar to the 
above expression with the fiuxes, probabilities and cross sections replaced by those for anti- 
muons. Similarly, the e~ event rate in a specific energy and zenith angle bin is expressed as 
follows: 



d^Ne 1 



dl^dE 27r 



dcosMEJ '^'^IdcosME^ '^'^ 



c"cc Dcff 



^Pq^ denotes the probability for transition from Vq, -^ vp. 



with the e+ event rate being expressed in terms of anti-neutrino fluxes, probabihties and cross 
sections. E and 9 in the above equations are true values of neutrino energy and zenith angle. 
We convert the above double differential event rates into those with respect to the measured 
neutrino energy Em and the measured zenith angle 9m- This is done using Gaussian resolution 
functions in energy and zenith angle, as explained in [10]. 

In the limit when only statistical errors are taken into account, the standard Gaussian 
definition of binned x^ is: 



Xstat 7 J 7 J ^ex V / 



Nf, 

i=Enibins j=cos6'mbins 'J 



Here, NfJ is the experimental and N^^ is theoretical number of events in the ij**^ bin. 

However, in addition to the statistical uncertainties, one also needs to take into account 
various theoretical and systematic uncertainties. In particular, our analysis includes a flux 
normalization error of 20%, a tilt factor [32] which takes into account the deviation of the 
atmospheric fluxes from a power law, a zenith angle dependence uncertainty of 5 %, an overall 
cross section uncertainty of 10 %, and an overall systematic uncertainty of 5 % [10]. These 
uncertainties are included using the method of pulls described in [32-34]. 

In this method, the uncertainty in fluxes and cross sections and the systematic uncertain- 
ties are taken into account by allowing these inputs to deviate from their standard values in 
the computation of Ny*^. Let the k**^ input deviate from its standard value by o"k ^k, where o"k 
is its uncertainty. Then the value of Njj*^ with the changed inputs is given by 

npuU 

N? = Nff (std) + £ cy ^k (10) 

k=l 

where N[|'(std) is the theoretical rate for bin ij, calculated with the standard values of the 
inputs and npull is the number of sources of uncertainty, which in our is case is 5. The ^k's 
are called the "pull" variables and they determine the number of cr's by which the k**^ input 
deviates from its standard value. In Eq. fllOp . cj' is the change in Njj^ when the k**^ input is 
changed by cTk (^-c by 1 standard deviation). The uncertainties in the inputs are not very 
large. Therefore, in Eq. (ITOl) we only considered the changes in Ny^ which are linear in ,^k- 
Thus we have a modified x^ defined by 

,, , ^ L^[f(std) + EkS^cHek-N^]' ^ 
X (Ck) = E W^ + E ^k (11) 

ij 'J k=l 

where the additional term ^^ is the penalty imposed for moving k*'^ input away from its 
standard value by ak ^k- The x^ with pulls, which includes the effects of all theoretical and 
systematic uncertainties, is obtained by minimizing x^('Ck), given in Eq. ( ITT]) , with respect to 
all the pulls .^k: 

xU = Mm^^ U'(^k)] (12) 



However, in general, the values of the mass-squared differences |Ani|^| and Am^i and 
the mixing angles 612, O23 and ^13 can vary over a range that reflects the uncertainty in 
our knowledge. To take into account the uncertainties in these parameters, we define the 
marginalized x^ foi' hierarchy sensitivity as. 



Y^ ■ = Min 

Amin 



x'(ek) + 



I ^^2^ I true _ |Am|i| 



a{\Aml\ 



sin^ 261'^^ - sin^ 2^13 



a (sin^ 29 



13 j 



(13) 



X^(^k) in the above equation, is computed according to the definition given in Eq. flTTl) . The 
following values for the errors are chosen 

• a (sin^ 2^13) = 0.02, 

• a (|Am2j) = 10% of |Am|i| = 0.25 x lO'^ eV^ 

In computing Xmin? ^^ varied the values of the neutrino parameters from their true values 
over ±3(7 ranges. Note that we have not done the marginalization over the mixing angle 6*23. 
In this work, we have considered only one value of ^23"° = tt/4. For this value of ^'23"'^, it is 
found that the x^jj^ always occurs at the true value itself. We do not marginalize over the 
solar parameters 612 and Am'^i in our computation, since they are known with very good 
precision [35]. 

We have computed the values of the x^ sensitivity to the mass hierarchy, choosing the true 
hierarchy to be normal. The true values of the neutrino parameters are taken to be 6'23 = 45°, 
e^2 = 33.8°, lAm^^l = 2.5 x 10"^ eV^, and Aml^ = 8.0 x 10"^ eV^ from the current 3a allowed 
ranges [35-38]. 

4 Results and Discussion 

We give our results for an exposure of 1 Mt yr (unless otherwise stated) for two input values 
of 6*13 and for detectors both with and without charge identification capability. When the 
detector has no charge identification capacity the x^ is defined as, 

{xl^)noID = Xl+fl (14) 

and 

{xll)noID = Xl+e- (15) 

That is, we sum over the particle and antiparticle events and then compute the x^- The total 
Xtot is the sum of xli and xli- 

When the charge identification capability is there, the x^ values for muon events are 
computed as 

Xl = xl + Xl (16) 
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where the individual x^ for muon neutrino and muon antineutrino events are evaluated and 
then summed. We assumed that there is 100% charge identification for muons throughout the 
energy range 1-10 GeV. That, however, is not true for electrons. In fact, there is no charge 
identification capability for electrons in the energy range 5-10 GeV and only a partial charge 
identification capability for lower energy electrons. Therefore the ~x^ for electron events is 
defined as 

Xe« = ^(Xe + Xg)l-5GeV^ + (1 " '^) (Xe+e)l-5GeV^ + (Xe+Js-lOGeF- (17) 

Here, n is the fraction of electron events in the energy range 1-5 GeV which are assumed 
to have exact charge identification capability. As mentioned earlier, we expect a ~ 20% 
probability of exact charge identification in the energy range 1-5 GeV. Hence we take ra = 0.2 
in our computation unless otherwise mentioned. In other words, 20% of the electron events in 
this energy range are assumed to have 100% charge sensitivity, while for the remaining 80%, the 
electron neutrino and antineutrino events are assumed to be indistinguishable. Realistically 
one expects a non-zero charge identification efficiency for the remaining 80% electron events 
in this range, which would give a better sensitivity to the hierarchy. Hence we have given 
a conservative estimate in our computation of the hierarchy sensitivity of a Liquid Argon 
detector. We also give an estimate of how the sensitivity would vary if n is varied from to 1. 
The third term in the above expression is the x^ contribution from the electron events in the 
5 energy bins in the range 5-10 GeV, for which charge identification is not possible. Hence 
this x^ value is calculated using the sum of electron neutrino and antineutrino events in these 
energy bins. XeZ represents the sum of all x^ contributions from electron events. 

In Tabled] we present the values of Xm ^^'1 Xe« with fixed parameters for two different values 
of the parameter 6*13. Here, y^^^i = Xm + xli is the total fixed-parameter hierarchy sensitivity of 
a Liquid Argon detector with charge sensitivity. In each case, (^-)^'^y^"-^ denotes the sensitivity 
with only statistical errors (Eq. ([9])) and (^^^^^pn/z denotes the sensitivity with theoretical and 
systematic uncertainties incorporated as pulls (Eq. ( TT2l) ). The values of (Xiot)^"" ^.re plotted 
against the input values of sin^ 26'i3 in Figure [TJ 

Figure [2] shows the variation of (Xtot)^^^'' (fixed-parameter) as n varies from to 1. 

In Table [2] we give the total values of the x^ hierarchy sensitivity with pulls and priors with 
marginalization over the neutrino parameters ^13 and Amgj^ for two different values of 6'^™^. 
These values, denoted by {Xtot)mim ^^^ fo^ ^ Liquid Argon detector with charge sensitivity 
and with an exposure of 1 Mt yr, and represent the principal results of this study. The same 
values are plotted against the input values of sin^ 2^13 in Figure [3l For comparison. Table [3] 
gives the marginalized results for hierarchy sensitivity for a Liquid Argon detector without 
charge sensitivity. 

5 Summary and Conclusions 

From the results described in the previous section, we note the following features: 

• Table [1] and Figure [1] give the values of fixed-parameter x^ for muon events with charge 
identification capacity and electron events with partial charge identification capacity as 
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0.04 
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0.10 


31.7 


14.7 


46.4 


25.1 


9.3 


34.4 



Table 1: Values of fixed parameter x^ without and with pull for a Liquid Argon detector (1 
Mt yr) with charge sensitivity. The contributions from muon and electron events are also 
shown separately. 
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16.2 
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Table 2: Values of total marginalized x^ with pull and priors, for a Liquid Argon detector 
(1 Mt yr) with charge sensitivity. Also shown separately are the contributions of muon and 
electron events to the total x^- 
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Table 3: Values of total marginalized x^ with pull and priors, for a Liquid Argon detector (1 
Mt yr) without charge sensitivity. Also shown separately are the contributions of muon and 
electron events to the total x^. 
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Figure 1: Values of fixed parameter x^ with pull versus the input (true) value of sin^ 2^13 for 
a Liquid Argon detector (1 Mt yr). Shown is (Xtot)^""- 

discussed in the previous section. We see that the values of fixed-parameter {XtotY^'''' ^^^ 
uniformly high (> 5a for both input values of 613) for a Liquid Argon detector with 1 Mt 
yr exposure. Since earth matter effects are proportional to sin^ 2^13, the x^ is expected 
to vary linearly with sin^ 26*1 



713- 



Prom Figure[21 we observe that the dependence of the total x^ on the charge identification 
capability for electron events is not drastic, since the contribution from the electron 
events in this energy range represent only about one-fifth of the total hierarchy sensitivity 
of a Liquid Argon detector. 

Comparing the marginalized results with and without charge sensitivity in Table [2] and 
Table [3], we can see that x^ in general increases after incorporation of charge identi- 
fication power of the detector. This is expected since the hierarchy sensitivity comes 
mainly from the matter effect which is different for neutrinos and antineutrinos. Hence 
a detector with charge sensitivity is better capable of probing this. Note that the muon 
X^ shows a sharper rise with the inclusion of charge sensitivity than the electron x^, 
since charge identification with high precision is possible over the entire energy range in 
the case of muons. 



From our study, a > 4cr sensitivity to the neutrino mass hierarchy is predicted for a 




20 40 60 80 100 

Percentage of electron events (1-5 GeV) with exact charge ID 

Figure 2: Values of fixed parameter ~x^ witli pull versus the percentage of electron events in 
the range 1-5 GeV taken to be with 100 % charge ID efficiency for a Liquid Argon detector 
(1 Mt yr). Shown is (xL)^""- Here, sin^ 2^13 = 0.1. 



magnetized Liquid Argon detector with 1 Mt yr exposure for sin^ 26'i3 > 0.05 or 6*13 > 6°. 
We redid our calculation for smaller exposure of 333 kT yr, which can be achieved 
by a smaller Liquid Argon detector (30 kT, 11 years) with charge sensitivity. In this 
case, we obtain the marginalized total x^ = 4.9, which is > la hierarchy sensitivity, for 
sin^ 26'i3 = 0.04, and the marginalized total x^ = 9.8, which is > 3(T hierarchy sensitivity, 
for sin^ 2^13 = 0.1. 

• For a detector without charge sensitivity, a 2.5cr marginalized hierarchy sensitivity for 
sin^ 2^13 = 0.04 is possible with 1 Mt yr exposure. 
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